Discrete Time Quantum Walk Approach to State Transfer 



Pawel Kurzyiiski^'^'l^] and Antoni Wojcik^ 

Centre for Quantum Technologies, National University of Singapore, 3 Science Drive 2, 117543 Singapore, Singapore 
^Faculty of Physics, Adam Mickiewicz University, Umultowska 85, 61-614 Poznan, Poland. 

(Dated: January 15, 2013) 

We show that a quantum state transfer, previously studied as a continuous time process in net- 
works of interacting spins, can be achieved within the model of discrete time quantum walks with 
position dependent coin. We argue that due to additional degrees of freedom, discrete time quantum 
walks allow to observe effects which cannot be observed in the corresponding continuous time case. 
First, we study a discrete time version of the engineered coupling protocol due to Christandl et. 
al. [Phys. Rev. Lett. 92, 187902 (2004)] and then discuss the general idea of conversion between 
continuous time quantum walks and discrete time quantum walks. 



INTRODUCTION 

Lattice models appear in quantum theory in various 
contexts. Just like most theoretical models in physics, 
they were first used to describe the properties of nat- 
ural materials, however nowadays we are reaching the 
stage when we are no longer bounded by nature, since 
our ability to engineer our own new systems having de- 
sired physical properties is developing very rapidly. This 
allows testing of abstract theoretical lattice models in the 
laboratory with a help of new extraordinary systems like 
graphene [T], cold atoms in optical lattices [2 , photons 
in arrays of linear optical elements 3J , or compounds of 
thin layers of different materials known as heterostruc- 
tures [4J. The new field of quantum information, which 
strongly relies on the capability of engineering and ma- 
nipulating of quantum systems, has also rekindled in- 
terest in lattice models. In particular, it has been shown 
that lattice models offer the possibility of universal quan- 
tum computation (5, t6j . 

In this work we consider two related fields of research: 
quantum walks on graphs [3 [S] and quantum informa- 
tion transfer in spin systems [9l llOj . In particular we are 
interested in the discrete time quantum walk realization 
of a quantum state transfer on a spin chain with position 
dependent couplings TT]. It is known that special cases 
of spin lattice models with one magnetic excitation can be 
interpreted as continuous time quantum walks [7]. Also, 
it has been shown that the relation between the two types 
of quantum walks is not trivial [T^l US] , because discrete 
time quantum walks require an additional degree of free- 
dom known as a coin. We show for the first time how one 
can convert the position dependence of couplings into the 
position dependence of coins. We also argue that discrete 
time quantum walks provide a more general framework 
for the description of quantum diffusion on regular lat- 
tices, since due to the extra degree of freedom one gains 
more control over the walk than is allowed in a standard 
continuous time scenario. Moreover, the evolution of dis- 
crete time quantum walks is much easier to simulate on a 
classical computer. Our approach is based on the relation 



between discrete time quantum walks and the Dirac-like 
equation to which the underlying quantum walk is trans- 
formed in the continuous limit of both time and space 

BASIC CONCEPTS 
State transfer in spin networks 

Consider a network for which every vertex corresponds 
to a spin 1 /2. We say that there is an edge between vertex 
i and j if the corresponding spins interact. Most authors 
choose to consider an XX model of spin-spin interactions 
[HI HO] , we follow them and set the Hamiltonian of the 
network to be of the form 

^= E (1) 

{i,i}es 

where E denotes the set of edges and cr* is a Pauli X ma- 
trix for the i'th spin. The above Hamiltonian conserves 
the total spin number along Z axis cr* , therefore one 
may restrict studies to a subspace with a fixed number 
of excitations, where by excitation we mean a spin point- 
ing up along Z axis. In particular, we are interested in 
an evolution restricted to a subspace with no excitations, 
which is invariant under the evolution generated by the 
Hamiltonian, and to a subspace with only one excitation. 
In the second case, the single excitation can be considered 
as a single particle walking on an underlying network. 

Next, imagine that there are two marked spins, a 
sender and a receiver, and that the whole network is in 
a ground state, i.e. in a subspace with no excitations. 
Then, the sender prepares his spin in a superposition 
ck| t) + /3| i), so that the whole network is in the state 

(a|t)+/3U»s® U;...;)net®U),.. (2) 

The goal is to employ the evolution generated by the 
Hamiltonian to transfer the quantum state from s to r 
in finite time T 

|i),®|U---i)«et®(«|t>+/3U»r, (3) 
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which is equivalent to transporting the excitation from 
spin s to spin r. 

There are many protocols which allow for a perfect, or 
a nearly perfect state transfer over spin chains. The basic 
techniques are: engineered couplings, wave packet encod- 
ing, and active control (see [HI HO] ) ■ Here we concentrate 
on engineered couplings. In particular, we will consider 
two protocols, namely the protocol of Christandl et. al. 
[TT] and a weakly coupled spin protocol [T71 [TS] . 



Continuous time quantum walks 

The state in continuous time quantum walks (CTQW) 
[7] is described by the position of a particle on a graph. 
The continuous evolution, which is governed by the 
Schrodinger equation, is determined by the Hamiltonian 
H which is an n X n Hermitian matrix proportional to 
the adjacency matrix of the underlying graph G which 
has n vertices 



vertex 



Hi 



^0, ifr(*,j)ei?(G); 
~ 0, else. 



(4) 



where E{G), as before, denotes the set of edges of G. 
The probability that at time t the walker, whose state 
is \^{t)) = e~^^'^\ip{0)), is located at vertex x is Px{t) — 
\{x\'4'{t))\'^ , where is the orthonormal basis spanning 
the space of vertices. 

For CTQW an a graph all of whose edges have equal 
weight, the corresponding Hamiltonian is simply the 
Laplacian of the underlying graph L — D — A, which 
is a diagonal matrix D minus the graph adjacency ma- 
trix A (whose entries Aij are either zero or one if vertices 
i and j are disconnected or connected, respectively). The 
diagonal entries Djj = dj denote the degree of vertex j. 
As a result, the sum of all elements in each column and 
each row of the Laplacian is zero. When the Laplacian 
plays the role of the Hamiltonian, the diagonal entries 
are important only in the case of irregular graphs. For 
regular graphs they are all equal, and since one can freely 
add or subtract from the Hamiltonian multiples of iden- 
tity without changing the dynamics of the system, they 
can be neglected. In case of a general CTQW on a chain 
the only nonzero elements of the Hamiltonian are Hjj^i, 
Hj^j^i and Hjj. In case of a single excitation subspace of 
a spin network, the elements of the corresponding CTQW 
Hamiltonian are given exactly by the strengths of spin 
couplings Hij — J, 



Discrete time quantum walks 

In discrete time quantum walks (DTQW) |8j the state 
of the system \x, c) is described by the position of a walker 
on a graph x and by the state of an auxiliary system c 
which determines the direction the walker is going to take 




FIG. 1: Due to the fact that interference phenomenon oc- 
curs for the second nearest neighbors even vertices can be 
interpreted as edges. 



in the next step. This auxiliary system is often referred 
to as the coin, which for a walk on a one-dimensional 
graph like a chain or cycle is simply a two-level system 
c =■<— , — In such cases, one step of the evolution is given 
by a unitary operator U = SC, which is a product of the 
coin operator 

C|x,— >) — cos6'|a;, - sinf^lx, ^), (5) 
C\x,^) = sm9\x,~?') + cosd\x,^), (6) 

and the conditional shift operator 

^|a;,^> = \x + l,^), (7) 
^|a;,^> = \x-l,^). (8) 

Therefore one step is given by 

U\x,^) = cos6'|x + 1,^) - sin6'|a; - 1,^), (9) 
U\x,^) = sin0|a; + l,^)+cos6'|a;-l,^). (10) 

The parameter 6 denotes the coin flip rate. Equivalently, 
the unitary evolution operator of one step can be written 
as 



U 



(11) 



where p is the momentum operator and aj are Pauli ma- 
trices acting on the coin space. 

In DTQW on chains and even cycles, interference oc- 
curs only between positions separated by a distance of 
two (second nearest neighbors), therefore one can con- 
sider two independent walks: one on even vertices and 
another on odd vertices. Due to this fact even vertices 
can be treated as edges as in Fig. [ij which will be cru- 
cial for our work. Since we are interested in the walk on 
vertices, not on edges, in our case one step of evolution 
will be given by U"^ and the initial state will be always 
supported on the vertex space only, as in the Ref. jl2| . 
Moreover, throughout this work we will consider DTQW 
on finite chains, which can be simulated by walks on cy- 
cles with a reversing coin at one vertex, i.e. a coin for 
which coin flip rate = 7r/2 — see Eqs. (5j6l and Fig. 

m 
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FIG. 2: DTQW on a finite chain is eflfectively simulated by 
a walk on a cycle with one reversing edge. 



CONTINUOUS LIMIT OF QUANTUM WALKS 

The continuous limit of DTQW on a chain results in 
the one-dimensional Dirac equation [11^16) . The opera- 
tor ( 11 1 can be rewritten as 



U = U{t) = f.^Vo.vt^^ay^t^ 



(12) 



where v is the velocity and w is the angular frequency. 
The two parameters are chosen so that for a unit time 



step vi ~ 1 and ujt = 9, thus recovering Eq. (11 1. Taking 



infinitesimal time steps and applying the Trotter theorem 
one obtains 



lim U{t/ny 



(13) 



The term in the bracket above corresponds to the Hamil- 
tonian of the system, which has exactly the form of the 
one-dimensional Dirac Hamiltonian 



H = cpaz 



mc^cr,, 



(14) 



The correspondence between (13 1 and (14) is fully estab- 



lished when one identifies v with the speed of light and w 
with the mass of the particle. In particular, ujt = 9, how- 
ever we take i = 1, therefore the mass of the quantum 
walker corresponds to the coin flip rate 9. However, due 
to the discrete nature of time and the periodicity of coin 
operator e*'""^, one is unable to differentiate between 9 
and 9 + fc27r. On the other hand, in the limit of infinite 
mass the particle should be almost immobile and the re- 
lation between the Dirac equation and DTQW requires 
9 ^ ^, hence one often takes m = tan 6* (see [Hj). Note 
that for small angles tan^ « 9. It is also important to 
notice, that due to the structure of the conditional trans- 
lation operator the continuization of time is intrinsically 
combined with the continuization of space, i.e. either 
both space and time are continuous, or both are discrete. 

The Hamiltonian of CTQW on a chain with uniform 
couplings is the following: Hi^i^i = Hi^i^i = —J and 
Hi^i = 2 J, however for regular graphs diagonal elements 



can be neglected and put equal to zero. The action of 
the Hamiltonian on a state localized at x is 



H\x) = ~J {\x + 1) - 2\x) + \x- 1)) . 

The above is a discrete version of the Laplacian. 
continuization of the above leads to 



H 



dx'^ 



(15) 



The 



(16) 



Obviously, this Hamiltonian describes a free particle in 
one dimension with J = therefore the coupling con- 
stant J can be interpreted as an inverse of the walker's 
mass. 

In the continuous limit, the quantum walker can 
be considered as a relativistic (DTQW), or as a non- 
relativistic free particle (CTQW). Main parameters gov- 
erning the dynamics of the walk, namely coins and cou- 
plings, are related to the particle's mass. Strauch [T2] 
showed that continuous quantum walk can be obtained 
from discrete one in the limit 9 ^ ^. Heuristically, this 
limit corresponds to a very heavy particle, for which it 
is much harder to observe relativistic effects, therefore it 
corresponds to a relativistic to non-relativistic transition. 
Another important feature of quantum walks is that in 
both models space is discrete, but only for DTQW time is 
also discrete. Once again, a heuristic explanation of this 
fact can be given as follows. Let the two quantum walks 
be discrete versions of Dirac and Schrodinger free parti- 
cle, respectively. In non-relativistic model time and space 
are treated separately, therefore there is nothing strange 
in discretizing only one of them. However, in relativis- 
tic physics time and space are combined and one has to 
discretize the whole spacetime at once, since it would be 
problematic to define consistent Lorentz transformations 
on spacetime which is only partially discrete. 



PERFECT STATE TRANSFER IN DTQW 

Let us apply the analogy between mass, coins and cou- 
plings to show that a perfect state transfer from one po- 
sition to another can be realized within DTQW. We con- 
centrate on the acclaimed perfect state transfer protocol 
introduced by Christandl et. al. [TT] . This protocol relies 
on properly engineered couplings, i.e. properly chosen 
terms of the Hamiltonian governing the corresponding 
CTQW. These couplings depend on position, therefore 
by our analogy, the corresponding particle has position 
dependent mass. On the other hand, mass is related to 
coin flip rate and so the corresponding DTQW should 
have position dependent coin. Below, we examine how 
the choice of such coins affects the dynamics of DTQW. 

The chain coupling between nodes n and n + 1 in 
the protocol is given by J„ = ^^Jn{N — n), where 
n = 1,2, . . . , N — 1 and A is a real parameter common 
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for all couplings. As a result, the effective mass of the 
particle for the transfer between positions n and n + 1 is 



1 



1 



(17) 



2J„ Xy/n{N^n) 
Next, consider a step operator of the form 

U = Uin)^e'P'''e"'y'^-, (18) 

where 0„ depends on position. The matrix form of the 
coin operator is given by 



Cn 



cos f, 
— sin^ 



sin 9n 

cos 9ri 



(19) 



Recall that we consider a double step operator C7^, since 
our walk on a chain of length N/2 is simulated by a walk 
on an even N-cycle for which even positions correspond 
to edges and we are interested in a walk on positions 
corresponding to vertices. Coins associated with edge 
positions correspond to mass, but there are also coins 
associated with vertices — odd positions of the cycle. 
Here, we set all of them to be identities (6*2^-1 = 0), 
although we note that a different choice (equal for all 
vertices) leads to similar results. The edge coins (even 
positions) are chosen according to the analogy between 
couplings and mass (to — tan 6) 



J 2k 



arctan ■ 



1 



2\. k 



' N_ 

K 1 



k) 



(20) 



where k = l,2,.. .,N/2 - 1. 

Our studies show that the transfer strongly depends 
on parameter A. For small values of A (large mass limit) 
the transfer is perfect. Around A = it starts to drop 
down and it recovers back to perfect transfer for A of the 
order 0(1) and greater (small mass limit). We show that 
the nature of transport for the two limiting cases is dras- 
tically different. In the small mass limit perfect transfer 
occurs due to dispersion-less nature of the operator e'*"^^ 
whose spectrum is linear, whereas in the large mass limit 
the generator of the operator is effectively given by a 
Hermitian operator whose action is the same as the action 
of the Hamiltonian introduced by Christandl et. al. |llj . 
Interestingly, in both cases we found that the transfer 
fidelity does not depend on the initial coin state. More- 
over, we checked numerically that a number of initial coin 
states, among which were eigenvectors of the three Pauli 
matrices, are always perfectly transfered. This fact al- 
lows us to conjecture that in DTQW case not only parti- 
cle, but also its intrinsic coin state is perfectly transfered: 

a|l,^)+/3|l,0 ^ a\N-l,^)+m-h^)- 

In Figure [3] we present the transfer of particle from 
position 1 to position TV — 1 for the corresponding DTQW 
with respect to evolution operator U'^. In this particular 
case = 30 and A = 0.03. As one can see transfer is 
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FIG. 3; Perfect state transfer from position 1 to position 
N — 1 for DTQW with position dependent coin (A*' — 30 and 
A — 0.03). Thin blue line — probability of occupying position 
1; thick purple line — pobability of occupying position A'^ — 1; 
dashed line — probability of occupying remaining vertices of 
the chain. 
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FIG. 4: Arguments (p of eigenvalues e^"^ of a position de- 
pendent coin DTQW double step operator for TV = 30 
and A = 0.03. Each dot corresponds to a doubly degenerated 
eigenvalue. 



perfect, and moreover the dynamics is periodic because 
after twice the transfer time the system goes back to the 
initial state, as in the continuous time case. 

The main idea of the authors of the protocol [Tl] was 
to find a mirror symmetric Hamiltonian with harmonic 
spectrum, which would provide periodicity of evolution 
and a perfect transfer between mirror symmetric posi- 
tions for times equal to half of a period. Numerical sim- 
ulations confirm that the spectrum of [/^ is also har- 
monic, yet it posses an additional relativistic property 
(see Fig. [4|. The arguments of eigenvalues of can 
be interpreted as quasi-energies. They split into two 
bands, which can be related to positive and negative en- 
ergy bands just like in the Dirac equation. Moreover, the 
two bands are separated by a band gap corresponding to 
the rest energy 2mc^. Due to the discrete nature of time 
for quasi-energies the band gap occurs twice, the first 
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centered around zero and the second centered around tt. 
Numerical simulations show that the width of the band 
gap is related to the inverse of A: in the limit A the 
band gap goes to tt, and for A — >■ oo the band gap goes to 
zero. This is as expect ed, since the corresponding mass 
see Eq. Ell. 



scales as ^ 



Large mass limit 



Recall the formulae for position dependent coin ( 19 ) 



and coin flip rate (20 1 and recall that arctan^ can be 
expressed as 



arctan^ — arcsin 



arccos 



1 



The entries of the coin matrix can be rewritten as 
2A7fe . „ 1 



cos f 2fc = 



v/rT4A^' 



sm U2k 



(21) 
, (22) 



where jk = y k — k) . In the large mass limit A — 0, 
which allows us to keep only the highest terms of Taylor 
expansion 



cos6'2fc « 2A7fe, sin 6*2/0 ~1- 



(23) 



The above limit resembles the one of Strauch [12 62k ^ 
|. The action of the double step operator [/^ on po- 
sitions corresponding to vertices (odd positions) can be 
approximated by 



C/2|2fc- 1, 
C/2|2fc- l,■ 



2A7fe|2fc+l,^ 
2A7fc_i|2fc-3, 



|2fc- 1,^), (24) 
+ |2fc- 1,^(25) 



Let us also introduce the following almost normalized 
states 



unitary operator can be approximated by [7 w 1 — iHt, 
where H is the Hamiltonian, therefore 

C/2|V'±(fc)) «z(l-ii/)|0±(fc)). (29) 

The above Hamiltonian is exactly the Hamiltonian of 
Christandl et. al., therefore the CQTW behavior is re- 
stored. Moreover, the action of H does not mix ± states, 
which is the reason why we observe not only the transfer 
of the walker, but also the transfer of its intrinsic coin 
state. 



Small mass limit 



Let us consider A — > 00. In this case a dispersion- 



less movement dominates the evolution, since Eq. ( 22 ) is 
approximately 



cos6'2fc«l, sin6'2fewO, 



(30) 



therefore the double step operator simplifies to U'^ ~ 
g-2ip(Tj ^ except at position N, where we apply the revers- 
ing coin which also causes a tt phase shift of the particle 
coming from the left. The transfer occurs in N/2 dou- 
ble steps. The initial state q;|1,->) + /3|1,^) evolves in 
the following way. The first part of this state, which is 
the right moving part, moves to the right, is refiected at 
position N and goes back to position — 1 to finish in 
state — q;|A^ — 1,^). The second left moving part goes 
left to = mod N, gets reflected and moves right to 
position — 1 to finish in state P\N — 1,^). This also 
works for any initial state localized at x, for which the 
destination point is N — x. In order to restore the final 
coin state to the initial state one has to apply the ay 
operation. 



Transition between the two types of transfer 



|V'±(fc)) = (|2fc - 1, ^) ± i\2k - 1, ^) 
±i2A7fe|2fc + l,^) -2A7fc_i|2fc-3,^)), (26) 

and normalized states 
1 



l</'±(fc)> 



|2A:-l,^)±i|2fc-l,^)) 



(27) 



Note that |0±(fc)) and \ip±{k)) are almost equal up to 
the factor of 0(A). The action of double step operator 
on \ip±{k)) yields 

U^\ij±ik)) = ii\^±ik))-i2Xjk\(t>±ik + l)) 

- z2A7fc_i|0±(fc-l))) + O(A2) (28) 

Now, it is enough to recognize that 2A7fc = Jk,k+i from 
the protocol [TT] and that for short times the action of the 



The unitary evolution operator U is the product of two 
operators S and C. Transfer properties rely on the inter- 
play between this pair. In the small mass limit C does 
not play any significant role, since it is close to identity. 
On the other hand, in the large mass limit C is crucial, 
because it makes the evolution similar to the one of the 
protocol [llj . As we pointed out in the beginning of this 
section, we observed that the transition between the two 
types of behavior occurs around A = Here, we give 
the heuristic explanation for this value. 

The time of transfer over a chain of length N for the 
protocol [TT] is T = J. In our case we consider chains 
which are effectively of length N/2, that is why we have 
to multiply A by two, which is evident in Eq. (20). As 
a result, the corresponding time of transfer should be 
given by T = ^ and indeed this is approximately what 
we observe. Next, let us recall that DTQW is related to 



relativistic quantum mechanics via its similarity to the 
Dirac equation. In relativistic physics a distance trav- 
eled by a particle in time t can be at most ct. Simi- 
larly, in DTQW in one step the particle can move only 
to neighboring positions, that is why in order to travel 
the distance from one end of the chain to another the 
particle needs to take at least N/2 steps. This gives the 
critical value of Xcr = above which the particle would 
move with velocity greater than the speed of light. Since 
this is impossible, a relativistic behavior has to domi- 
nate the evolution and as a consequence the transition 
between the two types of behavior can be interpreted as 
a non-relativistic to relativistic transition. 



Other techniques for perfect transfer 

As mentioned earlier, nearly perfect transfer in spin 
networks can be obtained in many ways. In principle, 
all transfer techniques can be employed in DTQW. For 
example, wave packed encoding (see [S] and references 
therein) is straightforward since one can prepare the same 
spatial wave packet in DTQW as the one considered in 
CTQW. However, position dependent coupling method 
is the one for which the correspondence between DTQW 
and CTQW is the most nontrivial. Bellow we would 
like to mention another example based on our coin-mass- 
coupling analogy. We present DTQW version of a pro- 
tocol in which two spins are weakly coupled to ends of a 
spin chain with uniform couplings |T7l [TB] . 

The CTQW Hamiltonian of a chain with weakly cou- 
pled ends is determined by coupling constants Ji,i+i = J 
for i = 2, . . . , — 2 and Ji^2 — Jn~i,n — aJ, where 
a <C 1. Due to weak coupling to the rest of the chain an 
interaction of the first and the last spin with the middle 
part can be treated as a perturbation. Without interac- 
tion the two spins would be in a degenerate state. Weak 
coupling breaks the degeneracy and perturbation theory 
predicts that there are new eigenstates which are even su- 
perpositions of two states: a state in which the first spin 
is up and one in which the last spin is up, or three states: 
the first spin up, the last spin up, or the middle of the 
chain being in one of its unperturbed eigenstates. The 
second case happens when the degenerated eigenvalue of 
the two spins is the same as one of the eigenvalues of 
the chain. Nearly perfect transfer is possible because ini- 
tial and final states are almost completely supported on 
these two/three states only (for the detailed discussion 
see [nillH])- 

In DTQW the above protocol can be realized in the 
following way. As before, we can simulate A^/2-chain on 
A^-cycle with one reversing edge. All coins are the same 
9k = 0, except 02 = 6*^-2 = f — e. We observe that 
DTQW behaves exactly like the corresponding CTQW. 
Moreover, we also observe the transfer of the coin state. 
In Fig. [5] we show the eigenstate population for A^ = 30, 
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FIG. 5; Eigenstate population for A'^ = 30, 9 — 7r/4 and 
initial state 11,—^). Left: Oi = 9n-i — 6 — tt/A. Right: 
62 ~ = 7i"/2 — £, where e — the initial state is 

almost entirely supported on four eigenstates only. 



p 




Time 



FIG. 6: Probability of population of positions 1 (solid) and 
A - 1 (dashed) in time. A = 30, e = f and e = . 

6 = J and initial state In case e 1 the ini- 

tial state is almost entirely supported on four eigenstates 
(two positive energy eigenstates and two negative energy 
eigenstates). In Fig. |6]we present the population of po- 
sitions 1 and A^ — 1 in time. The time of transfer scales 
asO(j^). 

We also observe another interesting phenomenon. If 
one changes the coin operator at position 2 from e*"^" 2 
to e'°'^ , the transfer is suppressed and the state is es- 
sentially localized at position 1 for all time. The reason 
for this behavior is that the initial energy at position 1 
differs from the energy the particle would have at posi- 
tion A^ — I, therefore the transport is forbidden due to 
energy conservation. In the next section we will explain 
this in further detail. 



HOW TO CONVERT BETWEEN CTQW AND 
DTQW? 

The conversion of CTQW into DTQW was considered 
by Childs |I3| . however his approach was purely mathe- 
matical and did not take into account the physical prop- 
erties of the underlying quantum walks. On the other 
hand, our approach is based on physical aspects of quan- 
tum walks, and hence by studying the DTQW version 
of a continuous process one can learn many facts about 
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the physics governing the corresponding system. In this 
section we show how to, in general, convert an arbitrary 
CTQW on a chain into DTQW. Once again, we are go- 
ing to implement walk on a A'^/2-chain using A'^-cycle with 
vertex and edge positions. 

Recall that in general CTQW on a chain the only 
nonzero elements of the Hamiltonian are iJjj+i, Hjj^i 
and Hjj. Diagonal terms, which are necessarily real, cor- 
respond to the potential energy 



u{j) = 77,- 



(31) 



Potential energy has been already considered in DTQW, 
mostly in the context of localization (see for example 
[T^). The kinetic energy for the transition between po- 
sitions j and j + 1, which is proportional to the second 
power of momentum divided by the mass, corresponds to 
the real part of the off-diagonal terms and Hj+ij. 

In this case Re[Hjj+i] = ^ . It was studied in the 



previous sections. The mass can be expressed as 



2Re[Hjj^i] 



(32) 



Finally, let us consider the imaginary part of Hjj^i. 
First of all, note that operator proportional to + 
1| + *|j + which gives rise to nonzero imaginary 

part of -ffjj+i, can be approximated as a first power of 
momentum operator for positions j and j + 1 (see [20j ) . 
Indeed, 

Y.-t\j){j + l\+Aj + l){3\ 
^2|j)(j|-|j+2)(j|-b-)(j- + 2|. (33) 

3 

The above operator resembles the Laplacian, which is 
the kinetic energy operator proportional to the square of 
momentum. The first power of momentum can appear 
in the Hamiltonian if there is a vector potential A acting 
on the particle. In this case we have to substitute (p)^ — > 

{v~ f 1 where e is the charge of the particle. As 
before, in the following, we assume c = 1 and we also set 
e — \. Using this analogy one arrives at 



ImfiJ. 



A 



(34) 



where in general both, the vector potential and the mass, 
are position dependent. The vector potential alone can 
be estimated as follows 



A 



1 

2 Re[H. 



(35) 



Our approach is based on the analogy between DTQW 
and the Dirac equation, whose general time-independent 



one-dimensional form is 

-i-^ - <^zij + m{x)ay%p = - u{x)^ tp. 

(36) 

At this point we can go to DTQW via application of the 
Trotter formula. One step of the corresponding DTQW 
is given by 

If — SCV — e^^'^' e^^"^^^^'^^^^^^^'^^^ e^'^^^\ (37) 

where j denotes position. To include the account of 
the scalar potential, we introduced a new operator V = 
Moreover, the new coin operator C contains the 



vector potential term. It is worth noting that vector po- 
tentials allow us to implement Hadamard coins using rel- 
ativistic analogy. Relativistic limit of Hadamard DTQW 
was considered in |16j. where the mass term in the cor- 
responding Dirac equation violated Lorentz covariance. 
Before we plug Eqs. pTl), (Isl and (psl) into Eq. (l37|. 



let us emphasize that in some cases it might be conve- 
nient to consider inverse tangent of both potentials, like 
we did previously with the mass, because of the period- 
icity of the arguments in the exponents. Moreover, since 
we consider division of positions into edges and vertices, 
scalar potential terms should act on vertices, whereas 
mass and vector potential terms should act on edges. As 
a result the double step operator can be taken as 



[/2 — SCSV — e*P'^-'e*(™^"'^'*^'^'' 



lCi)g*I«' 



(38) 



where A 



arctan I 



arctan 
and 



( 1 Im[gj., + i] 

— arctan(_ff. 



2Re[J?j,j + i] 

A similar conversion can be done for higher- 
dimensional quantum walks. Moreover, it is clear that 
using the above analogy one can convert the correspond- 
ing DTQW back into CTQW. However, there are certain 
kinds of DTQW which cannot be easily converted into 
CTQW, therefore let us now concentrate on the reverse 
problem of conversion of DTQW into CTQW. An exam- 
ple of a case in which the conversion is problematic is 
a DTQW with position dependent coin, where not only 
the coin flip ratio 6j depends on position, but also the 
generating operator Uj changes from position to position. 
In the previous section we considered DTQW version of 
a weakly coupled spin protocol and observed that the 
transfer can be suppressed due to a change of the first 
and the last coin flip generators from ay to ■ The coin 
degree of freedom is a qubit, therefore the coin generator 
is a linear combination of (Jy and CTz. In the Bloch 
sphere picture, coin operation is a rotation of a qubit 
about an axis n, where the generator of rotation is n • (?. 
In general, n can point in any direction on a surface of 
a three-dimensional sphere, however in the case of the 
Dirac equation (36) it is bounded to the yz-plane. Due 



to this fact, simple conversion from DTQW into CTQW 
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is not possible. As a result, DTQW offers a greater pos- 
sibility of control over the transport than is allowed in 
the CTQW case. 

Up to now we have considered only one-dimensional 
graphs, like chains. In two, and higher, dimensions one 
encounters even more complicated obstacles. Firstly, the 
typical formulation of a DTQW on a d-dimensional grid 
uses a 2(i-dimensional coin. Each coin dimension cor- 
responds to a different direction for the walker to take 
in the next step. This is in contrast to the Dirac for- 
mulation of relativistic motion, where for example in 
the three-dimensional case the coin space is only four- 
dimensional. The conditional shift operator of the corre- 
sponding DTQW is of the form 



Exp i ^ pj (g) 



1] 



(39) 



where j = 1, 2, . . . , d enumerates different orthogonal di- 
rections, 7j = |2j)(2j| — |2j — l)(2j — 1| acts on the coin 
space and pj is the corresponding momentum operator 
for direction j. Note that [7j,7fc] = for all values of 
j and A:, whereas the corresponding Dirac matrices anti- 
commute {7j,7fe} = 2^^^. Moreover, 7| = 1 for all j, 
whereas for DTQW J^jlj — ^- Due to this reason, 
one is unable to recover the relativistic energy relation 
E"^ = \p\'^+m'^, but rather obtains E'^ = J2j Ictjl'^p'j + ^ 
where J^j I'^iP — 1 ^-i^d is a parameter related to 
properties of the relevant coin operator, which can be in- 
terpreted as the mass. Since our method of conversion 
is closely related to the similarity of quantum walks to 
the Schrodinger and the Dirac equations, it is somehow 
expected that problems with conversion of DTQW into 
CTQW appear when the description of the DTQW de- 
parts from the Dirac one. 

It is important to mention here the research on quan- 
tum walk search protocols on d-dimensional grids. The 
efficiencies of the two protocols, the DTQW pTJ [22] and 
the CTQW [U, are quite similar. Only in dimensions 
d = 2,3 the DTQW protocol is faster than the corre- 
sponding CTQW protocol. In [20], the authors of the 
CTQW protocol improved the efficiency of their model, 
however this was done at the cost of an additional coin 
degree of freedom, which had to be included — departure 
from the standard CTQW structure was inevitable. The 
heart of the DTQW protocol [5T] is the so called Grover 
coin, which for d = 2 is of the following form 



1 1 1 



1\ 



-1 1 1 
1 -1 1 
1 1 -1/ 



(40) 



For the two-dimensional Grover walk the probability dis- 
tribution remains localized at the initial position for all 
time. This effect was explained by Inui et. al. |24| . 



who noticed that eigenvalues of the Grover walk are 
highly degenerated. More than half of the eigenvalues 
are ±1. There is no CTQW analogy of the Grover 
walk. Since the Grover walk obeys translational symme- 
try, the corresponding CTQW version should also pos- 
sess this property. However, for a translational sym- 
metric CTQW on a iV X TV grid with periodic bound- 
ary conditions, coupling constant J = \J\e^^ and di- 
agonal term A, the set of eigenvalues is given by 
A + 2\J\ (cos(<^ - 2TTk^/N) + cos((/? - 2'Kky/N)), where 
kx,ky = 0,1, . . . , N — 1, and it is clear that there is no 
possibility of such a high degeneracy. 



CONCLUSIONS 

In this work we studied the DTQW version of the 
quantum state transfer, which is originally described as 
a continuous time process. While in the continuous time 
scenario the transfer properties depend on couplings be- 
tween neighboring positions, in the discrete time case it 
is the coin operator which is responsible for the perfect 
transmission from one position to another. We applied 
the analogy between the DTQW and the Dirac equa- 
tion and between the CTQW and the Schrodinger equa- 
tion to show that both coins and couplings can be in- 
terpreted as the mass of a particle, which allowed us to 
transform tje CTQW into the DTQW. We examined the 
DTQW versions of perfect state transfer protocols stud- 
ied in [TTl HZl [IB]- We found that in DTQW versions 
not only the particle, but also its intrinsic coin state is 
perfectly transfered. The general method of the CTQW 
transformation into the DTQW has been also discussed. 
Finally, we argued that some DTQW's do not have the 
corresponding CTQW versions due to the fact that in 
some cases the dynamics of DTQW is much richer than 
the CTQW one. 

One of the biggest advantages of CTQW is that it can 
be easily defined on any graph, whereas DTQW is nat- 
ural on regular graphs only, where the same coin degree 
of freedom can be used for all vertices. However, many 
problems in physics, as well as in computer science, are 
defined on regular graphs. The above discussion shows 
that in many cases DTQW is more general than CTQW, 
since it allows to observe some effects which cannot be 
observed in CTQW. It would be interesting to find phe- 
nomena which can be observed in CTQW only. 
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